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Abstract 

We have explained in detail why the canonical partition function of Interacting Self Avoiding Walk 
(ISAW), is exactly equivalent to the configurational average of the weights associated with growth walks, 
such as the Interacting Growth Walk (IGW), if the average is taken over the entire genealogical tree of the 
walk. In this context, we have shown that it is not always possible to factor the the density of states out of 
the canonical partition function if the local growth rule is temperature-dependent. We have presented Monte 
Carlo results for IGWs on a diamond lattice in order to demonstrate that the actual set of IGW configura- 
tions available for study is temperature-dependent even though the weighted averages lead to the expected 
thermodynamic behavior of Interacting Self Avoiding Walk (ISAW). 

PACS numbers: 05.10.Ln, 36.20.-r 
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I. INTRODUCTION 



Monte Carlo simulations, based on Metropolis sampling algorithm, have immensely con- 
tributed to^ur understanding of a variety of complex physical systems and their thermodynamic 
behavior [1]. Yet, they are also known to be computationally inefficient in situations where low 
entropy microstates of a canonical system are to be sampled in sufficient numbers in order to ob- 
tain accurate estimates of relevant thermodynamic parameters. Employing suitable bias sing rules 
that favor these microstates, especially in irreversible growth models, does not ensure a neat solu- 
tion to this problem because the associated weights could wildly fluctuate leading to large errors. 
Wang-Landau flat histogram algorithm [|20 is a recently proposed dynamical method for efficiently 
sampling such microstates and hence, for accurately computing the Density of States (DoS) of the 
system under study. 

A particularly interesting idea, highlighted recently by Prellberg and Krawczyk [|3|] in the con- 
text of Self Avoiding Walk (SAW) M, is to recognize the DoS as simply the average of weights 
associated with all Monte Carlo attempts to sample the required microstates; acceptance of a mi- 
crostate depends on how the associated weight compares with the average and is implemented in 
such a way that the energy histogram becomes progressively flatter. In this step-by-step growth 
model, also known as the Kinetic Growth Walk (KGW) [jSD, the weight associated with a mi- 
crostate (a walk configuration) is the product of single step weights which are, in fact, the local 
microcanonical partition functions or equivalently the number of available directions for the in- 
dividual steps. Since the model is athermal, it is intuitively clear that the average weight could 
provide an estimate of the DoS. 

If, on the other hand, the individual steps are sampled on the basis of the energy being gained 
by the walk, their associated weights are customarily taken to be equal to the inverse of the cor- 
responding step probabilities; again the configurational average of their products could lead to an 
estimate of the DoS. An interesting alternative is to set the single step weights equal to the local 
canon/ca/ partition functions, a straightforward generalization of the athermal case. Assuming that 
the configurational average of their products over all Monte Carlo attempts leads to an estimate 
of the canonical partition function, it is not necessary that it would be a sum of terms that are 
factorizable into the (athermal) DoS and the corresponding Boltzmann factor. 

In this paper, we clarify this point by using Interacting Growth Walk (IGW) model [jg], which 
is a finite temperature generalization of the KGW. We explain how the canonical partition function 
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for walks of given length is exactly equal to the average of the products of local partition functions 
associated with all possible walk configurations including those that have failed to grow to the 
full length. In this case, we show that the density of states (DoS) cannot be factored out of the 

canonical partition function. Further, we present Monte Carlo results for IGWs on the diamond 
lattice and show that the set of IGW configurations available for study is temperature-dependent, 
in contrast to the athermal set of KGW configurations, even though the weighted averages lead to 
the expected thermodynamic behavior of SAW. 

II. STEP-BY-STEP GROWTH OF A SELF-AVOIDING WALK 

Consider a SAW configuration, Ck-i, made up of an ordered set of K — 1 directions, 
{fii, fi2, - ■ ■ , fJ-K-i}, taken consecutively on a regular lattice of coordination number z. Recog- 
nizing the fact that self-avoidance for the current step is a non-local requirement involving the 
entire walk configuration, we denote by ax (Ck-i) the number of acceptable or available direc- 
tions for the K*^ step. The walk proceeds further only if the direction chosen for the next step is 
acceptable; it is 'trapped' at the {K — 1)*'* step and will not grow further, if axiCK-i) — 0. In 
other words, local acceptability criteria decide the length and configuration of the walk. 

It is quite likely that an acceptable step leads to a site some of whose nearest neighbours are 
sites through which the walk configuration has already grown. Let n{iJ,K', Ck-i) be the number 
of such non-bonded nearest neighbours, also called contacts, encountered by the K^^ step in the 
direction jix- Clearly, the total number of contacts in the configuration Ck is the sum of contacts 
made by each step in the walk. By assigning a quantum of energy, say e, to each of these contacts, 
we will be able to treat the walk as a thermal object. 

There are z possible directions for the first step taken from an arbitrary lattice site, called the 
'origin'; each of these directions leads to {z — 1) possible directions for the second step, and so 
on until the first steps of the walk are taken without making any contact. The total number of 
possible configurations identified upto steps is therefore Bm^ = z{z — 

Some of these configurations would make contacts, for the first time, at the (M^-l- 1)*'* step. For 
example, 4, 2 and 1 on honecomb, square and triangular lattices respectively. Subsequently, 
each of these contact making configurations would have, say, aM^+2(CM^+i) directions for the 
[Mz + 2)*^ step. Some of these growing configurations will be geometrically 'trapped' at the Nf- 
step and hence caimot grow further. The minimum value of A^^ at which trapping occurs for the 
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first time depends on the lattice on which the walk is grown. For example, minAT^ — 9,7 and 6 on 
honeycomb, square and triangular lattices respectively. 

If we map out all the configurations till they are either geometrically trapped or have grown 
to their full specified length, say N, we have the genealogical tree, Z^, of SAWs of length less 
than or equal to N. Clearly, any trapped iT-step SAW configuration {K < N), realizable by a 
growth algorithm, is a branch of the genealogical tree. It flowcharts all possible outcomes of a 
growth algorithm and is, in fact assumed to be made up of all possible A^^-step SAWs as well as all 
possible trapped X-step SAWs {K < N). The shortest branch of the tree is of length N^, at which 
length the growth of the branch is stopped on a given lattice. Counting the number of branches of 
length K in the tree, Zjv, for all K in range [A^^, A^] could provide a complete description of the 
treelt is of interest to know whether a growth algorithm that uses local growth rules to trace out a 
branch can estimate, in particular, the number of branches of length A'^, which is equal to the total 
number of A^^-step SAWs, Z^. 

Let P{Ck, P) be the probability of growing a X-step walk, or equivalently of realizing a in- 
step branch of the genealogical tree Zn, at an inverse temperature p. Since any branch of this 
tree is a either a realizable SAW configuration of length equal to A?^ or a realizable trapped SAW 
configuration of length less than A^ , and vice versa, we have the normalization for their growth 
probabilities: 

f2 (j2piCK,P)]^i (1) 

K=N, \Ck / 

Since Zn (N > N^) does not have a branch of length less than A^^, P{Ck<n^, P) — ^ and hence 
the summation is from A^^ onwards. 

In terms of the single step probabilities, Cl-i,P), we can write 

Dfr ff^ J nf=M.+i Pil^L-, Cl-1, P) if all the steps are acceptable 

I if L*" step is not acceptable(L < K) 

where the prefactor B]^^^ = [z{z — is the probability of growing the initial segment of 

length < N^, which does not make any contact during its growth and hence is athermal. 

Single step probabilities, Cl-i, P), are in general assumed to be temperature dependent 
so as to take into account the possibility of bias due to contacts. However, in the case of KGW, it 
is temperature-independent and is given by 

p{lJiL]CL-i) ^ — r; ai(Ci_i)>0 (3) 



It is locally normalized over all acceptable steps. In the case of IGW, on the other hand, 
p{fj,L', Cl-1, (3) is a temperature-dependent locally normalized jump probability for the L*^ step 
in the direction /x^: 

g/3n(/i£;C£_i)e 

p{fiL]CL-i,f3) = — -; ai(CL„i)>0 (4) 

where ai(CL_i) is the number of available sites for the L*'^ step and, without loss of generality, e 
may be set equal to unity. 

III. ESTIMATING THE CANONICAL PARTITION FUNCTION USING THE GROWTH PROB- 
ABILITIES 

Growing a SAW of length more than the minimum steps with history-dependent, local 
step-probabilities has an important consequence - namely, that its growth probability, P, will not 
be the same as that obtained by growing the same configuration in the reverse order. The reason 
for this is the different local growth environments in which the walk has to sample its next step. 
One way of quantifying and correcting for this is to assign a canonical weight to the L*'* step in 
the direction fi^ of a growing configuration Ci_i, if acceptable {i.e., if Cl_i) = 1), 

/3n(^i;Ci_i) 

u;{f,L-CL-uf3)^-. ^ ^; 1<L<N (5) 

It may be recognized immediately as Grassberger's PERM-B weight OD for IGW and is, of course, 
zero for a step that is not acceptable. These weights assigned to the individual steps, in turn, define 
a canonical weight that can be assigned to an A^-step configuration, Cn'- 

N 

Wm{Cn,P) = Bm^ n (6) 

g/3n(Cjv) 

where the prefactor Bm-^ = — is due to the fact that the walk does not make any contacts 
for the first steps, and u^Cn) ( = J2l=Mz+i "^il^L'-, ) is the total number of contacts in an 
A^-step configuration C^. Note that Wn{Cn, (3) = 0, by definition, for walks of lengths less than 
N. 



5 



The average value of Wn{Cn, /3), taken over the entire genealogical tree Z^, is given by 

= J2Wn{Cn,P)P{Cn,P) (9) 

The last identity is due to the normalization of the growth probabilities given by Eq.(l), as well as 
due to W^iC^, P) = 0, by definition, for walks of length less than A^. From Eq.(7), it is clear that 

(W^iv(/5)) = 5^e^"("-)^Z^(/?) (10) 

where Zjy^) is the canonical Partition Function for A^-step Interacting Self Avoiding Walks 
(ISAW) [8]. In other words, averaging the canonical weights, defined in Eq.(7), over the entire 
genealogical tree, Zjq, gives the exact value of Zn{[3). 

Because mapping the entire genealogical tree, Zn, is quite a formidable task for large values 
of N, Monte Carlo methods are employed for estimating the average value of PV7v(/3): 



Wn^s = " ^ (11) 

where the summation is over all successful A^-step walks and S is the total number of attempts 
made to generate them. With the probability of growing a configuration implicitly taken care by 
the algorithmic rules, this equation may be recognized as an equivalent of Eq.(8); the Monte Carlo 
estimate of W^iP) will approach that given by Eq.(lO) as the number of attempts, S, becomes 



sufficiently large (i.e., Wn{(3)s^^ {WN{f3))). 

It must be noted that the identity, Eq.(lO), holds good whatever be the growth walk used. For 
example, the single-step probability of KGW is temperature-independent (Eq.(3)) whereas that of 
IGW is temperature-dependent (Eq.(4)), yet either of them could be used for estimating Zm{P). 
In order to understand the basic difference between them, we first rewrite Eq.(9) in the form, 

nx{N) I \ 

{WnW)) = Yl $^W^^(CjV,n,/9)P(CjV,n,/?) (12) 
n=0 yCjv,„ J 
nx(N) 

= Yl (^N,nm (13) 
n=0 

where the inner summation is over all the A^-step configurations that make n contacts and nx (N) is 
the maximum number of contacts an A^-step configuration will make. A comparison with Eq.(lO) 



immediately leads to the identification, 

nx(iV) nx{N) 

J2 (WnAP)) = Z^W) = J2 9N{n)e^'' (14) 

n=0 n=0 

where QNin) denotes the total number of configurations making n contacts, or equivalently the 
Density of States (DoS). This begs the question whether a term-by-term identification is also im- 
plied. 

In the case of KGW, the weights being temperature-independent, we have the identity, 

{WN,n) = Yl WN{CN,n)P{CN,n) = 9N{n) (15) 

CjV,n 

that leads to the factorizability of the individual terms (Wjv,n(/5))> namely, 

{W^^M) = (V^iv,n)e^" (16) 

In other words, term-by-term identification in Eq.(14) is meaningful for KGW because the ather- 
mally generated configurations are assigned appropriate Boltzmann factors a posteriori. 

On the contrary, the probability of generating an IGW configuration is temperature-dependent, 
and from Eqns. (4 - 7), we have the product. 



N /aL{CL-i) 



U{n-n{CN)) (17) 



where 5-function ensures that the configuration Cn has n contacts. It is clearly not factorizable 
into temperature-independent (DoS) and temperature-dependent terms, like in the case of KGW 
(Eqs.(15 and 16)). Yet, Eqs.(10 and 12) ensure that the configurational average of IGW weights 
(or equivalently, the PERM-B weights) provides an estimate of the canonical partition function 
and hence can also be used for computing the energy fluctuation as a function of temperature. 
A Monte Carlo estimate of {WN,n{P)) follows from Eq.(ll): 



W^,Ms = g ■ (18) 

In the case of KGW, the weights are temperature-independent and we immediately have a Monte 
Carlo estimate of the Density of States, QNin), from Eq.(15): 



9N{n)s = -^^^ (19) 



In the case of IGW, however, such a direct estimate of the DoS is not possible. 



A dynamical algorithm for estimating Wk{(3)s any given step K = 1,2, - ■ ■ , N has re- 
cently been proposed by Prellberg and Krawczyk ||3|]. Let Wk{Ck, P) be the canonical weight 
assigned to the {S + 1)*^ configuration, Ck, at the K^'^ step. This configuration is either allowed 
to grow further or is terminated depending on how its weight compares with the updated average, 



WK{f3)s+^ = [SWK{f3)s + Wk{Ck, f3)]/{S + 1): 

Wk(Ck,/3) I <1 continue the walk with probability r 

vk^Ck, P) = = < (20) 

^k{P>s+i I >1 continue (enrich) the walk 



Since WKiP)^ is an estimate of the canonical partition function, Zk, the parameter 
rK{CK,P) could be considered as the Boltzmann factor e^^^^'^^''^\ where AF(Ck,P) = 



(3 ^\og\WK{CK,l3)/ VFx(/3)g^-^]. Hence, without enrichment, Eq.(20) is a simple Metropolis 
criterion that is also physically meaningful for the growth walks. 

IV. A MONTE CARLO EXAMPLE 

Using this simple Metropolis criterion, without enrichment, we have estimated, and shown 
in Fig.l, normalized fluctuations in the number of contacts per monomer, a'^{m)/N, for fairly 
short IGWs of length upto 256 steps on a diamond lattice. As the walk length increases, the peak 
shifts towards the expected value, /5 ~ 0.44 [^, indicating thereby that the specific heat data 
are not sensitive to the factorizability of the canonical weights associated with KGW and IGW 
respectively. 

Large fluctuations in the weights could be a severe problem, as illustrated in the inset of Fig.l 
for example. In fact, sample loss due to pruning outweighs the reduction in loss due to attrition 
at higher values of [3\ consequently, reduction in the effective sample size may also contribute 
to wild fluctuation in the weights. The data presented in the inset correspond to IGWs grown at 
(3 = 0.7; it took about three hours of Pentium IV 3.2 GHz processor time to collect them. While 
the unweighted variance per monomer increases monotonically with the length of the walk, the 
weighted ones fluctuate wildly; nevertheless, a definite trend could still be discerned by a simple 
adjacent averaging of the weighted data. 

The overflow problems plaguing the IGW weights for long walks are taken care of by normaliz- 
ing the canonical partition function for the current step with the effective coordination number at a 
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FIG. 1: Normalized weighted variance, a'^{m)/N, per monomer as a function of the inverse temperature, 
/3 for IGWs on the diamond lattice for lengths N = 64, 128 and 256. The peak shifts towards the expected 
value f30 ~ 0.44 for longer walks. INSET: Normalized variance per monomer plotted as a function of the 
length of the walk grown at /? = 0.4. The unweighted data is monotonically increasing with N, whereas 
the weighted ones are highly fluctuating. A simple 10-point adjacent averaging of the unweighted data is 
shown by the continuous line. 

given value of P, which is nothing but the asymptotic estimate of the configuration-averaged, local 
canonical partition function per step. It is observed that the weights are log-normally distributed 
(Fig. 2), and its peak value, {ln[W{N, P)]), is asymptotically proportional to for any given value 
of p. In fact, we find that {ln[W{N,p)]) = a^(/5)+6^(/?)yiV+c^(/?)A^ (upper left inset of Fig.2), 
where the small scale-shift parameter, a^(/5), may be ignored for large N; /3-dependence of the 
parameters and Cw is presented in the right inset of Fig. 2. So, in the asymptotic limit (N oo), 
{ln[W{N, P)]) (X N which is consistent with the expected scaling form {W{N, P)) ~ z^jj, where 
Zeff is the effective coordination number of the lattice. It must, of course, be mentioned here that 
(/n[Vr(A^, P)]) 7^ ln[{W{N, P))], even though they both vary linearly with N for long walks, as 
demonstrated in the lower left inset of Fig. 2 for P = 0.5. 

On the other hand, the weighted distribution for the number of contacts per monomer, P{m/N), 
is quite sensitive to the the type of growth walk used. This is evident from Fig. 3 in which the 
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FIG. 2: Normalized distribution of logarithm values of the IGW (or PERM-B) weights for N = 256 at 
/5 = and 0.8. UPPER LEFT INSET: Peak values of the distribution,(/n[H^(Af, /?)]), as a function of //V 
for P = 0,0.4 and 1.0 (bottom to top). The data fall on a parabola, {ln[W{N,P)]) ~ b^VN + c^iV. 
RIGHT INSET: The parameters, 6^ and , versus /3. LOWER LEFT INSET: A demonstration of the fact 
that {ln[W{N, (3)]) / ln[{W{N, (3))] at /5 = 0.5. The deviation becomes larger at higher values of N. 

weighted P{m/N) for IGW grown dX (5 = 0.4 is compared with the distribution for KGW canon- 
ically reweighted at /3 = 0.4. It is clear that the IGW configurations are more compact than 
the KGW configurations, eventhough they both have been counted as ISAW configurations at 
(3 = 0.4. The distribution has a non-Gaussian functional form, P{m/N) ~ —{m/N)°', espe- 
cially away from the peak region, as is evident from the inset of Fig. 3. This is in agreement with 
Baumgaertner's analysis 190. 

V. SUMMARY 

We have explained in detail how the configurational average of the weights associated with 
growth walks, such as IGW, provides an esimate of the canonical partition function. We have 
drawn attention to the fact that the individual terms of the canonical partition function are not 
expressible as products of the density of states and the Boltzmann factor if the local growth rule 



10 



0.06 



0.04 - 

z 

E 

oT 0.02 
0.00 



KGW to ISAW 






p = 0.4 
N = 256 



0.0 



0.1 



0.2 0.3 

m/N 



0.4 



0.5 



FIG. 3: Normalized probability distribution of the fraction of contacts for walks of length N = 256 at 
P = 0.4. The one peaking at a smaller value of m/N corresponds to KGW canonically reweighted to 
P = 0.4. The other one peaking at a higher value of m/N corresponds to IGW grown at /? = 0.4 and 
(PERM-B) weighted. INSET: Same data, but ln[—ln[P{m/N)]] is shown as a function of ln[m/N]. 

is temperature-dependent. In fact, the parameter, P, that tunes the step-by-step growth of an IGW 
may be interpreted a posteriori as the inverse 'bath' temperature if the growing configuration is 
assigned a PERM-B weight, as given by Eq.(5). Numerical support for this is provided by the 
specific heat data. Yet, the actual set of configurations generated depends very much on whether 
the growth rule employed is athermal or not. For example, at large values of P, the IGW growth 
rule may not generate certain (compact) configurations, which are expected to be realizable if P^^ 
were a bath temperature. We have demonstrated some of these subtle points by presenting Monte 
Carlo results obtained for IGWs on a diamond lattice. 
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